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Abstract 

Harmonic oscillator in Fock space is defined. Isospectral as well as 
polynomiality-of-eigenfunctions preserving, translation-invariant 
discretization of the harmonic oscillator is presented. Dilatation- 
invariant and polynomiality-of-eigenfunctions preserving discretiza- 
tion is also given. 
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Undoubtedly the harmonic oscihator plays a fundamental role in science. 
Goal of present note is to give a review of different representations of quan- 
tum harmonic oscillator and its deformations in terms of the elements of 
Heisenberg algebra, differential, finite-difference, discrete operators. 

The Hamiltonian of harmonic oscillator is defined by 

^ = -£+'"'^'' 
where oo is the oscillator frequency. The eigenfunctions and eigenvalues are 
given by 

^k{x) = //fc(V^x)e-'^'^ , Ek = u{2k + 1), A; = 0, 1, . . . (2) 

where Hk is the kth. Hermite polynomial in standard notation. Without 
a loss of generality all normalization constants we put equal to 1. The 
Hamiltonian (|^) is Z2-invariant, x — x, which leads to two families of 
eigenstates: even and odd, symmetric and anti-symmetric with respect of 
reflection, correspondingly. This property is coded in parity of the Hermite 
polynomials: 

H2n+p{^x) = xPL^r~^\oJx^) , n = 0, 1, . . . (3) 

where L^n'\y) is the nth associated Laguerre polynomial in standard nota- 
tion, and p = 0,1 has a meaning of parity. Hereafter we can call a ground 
state the lowest energy state of parity p: 

= xPe-'^^ , (4) 

Thus the formula (4) makes an unification of both possible values of parity 
and for the sake of simplicity we will call (4) the ground state eigenfunctions 
without specifying parity. 

Make a gauge rotation of the Hamiltonian (1) taking a gauge factor the 
ground state eigenfunction (4) and change variable x to y = iox"^, which 
incorporate the reflection symmetry. Finally, after dropping off the constant 
terms we get an operator 

h{y,dy) = -- {^l^\x))-'n^l^\x)\y^^,2 = Aydl-A{y-p-\)dy 

^ (5) 

with the spectrum (— 4n), where n = 0, 1, 2 ... . The operator (5) simulta- 
neously describes a family of eigenstates of positive parity if p = and a 
family of eigenstates of negative parity \ip = 1. We will call (5) the algebraic 
form of the Hamiltonian of the harmonic oscillator. The word 'algebraic' 
reflects the fact that the operator (5) has a form of linear differential oper- 
ator with polynomial coefficients and furthermore possesses infinitely-many 
polynomial eigenfunctions. The latter implies that any eigenfunction can be 
found by algebraic means by solving a system of linear algebraic equations 
The algebraic form (5) admits a generalization of the original Hamilton- 
ian (1) we started with. If we assume that the parameter p can take any 
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real value, p > —1/2, one can make a inverse gauge transformation of the 
operator (5) back to the Hamiltonian form and we arrive at 



4yd^^-4{y-p-^)dy 



y-p/2gy/2 



-Uk , (6) 



which is known in literature as Kratzer Hamiltonian. It is worth to mention 
that this Hamiltonian coincides also with 2-body Calogero Hamiltonian. 
Hereafter we will call the system characterized by the Hamiltonian (6) the 
harmonic oscillator. 

The resulting Hamiltonian (6) is characterized by the eigenfunctions 
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^fc(x) = xPL„^"^(u;x^)e-"^ , (7) 

which coincides with (2) at p = 0, 1. The spectrum (6) is still equidistant 
with energy gap lo and after appropriate shift of the reference point it co- 
incides with the spectrum of the original harmonic oscillator (1). Thus, the 
deformation of (1) to (6) is isospectral, which is, of course, well-known. 

In order to move ahead let us introduce a notion of the Fock space. Take 
two operators a and b obeying the commutation relation 

[a, b] = ab-ba = I, (8) 

with the identity operator / on the r.h.s. - they span a three-dimensional Lie 
algebra which is called the Heisenberg algebra h^, ■ By definition the universal 
enveloping algebra of /13 is the algebra of all normal-ordered polynomials in 
a, b: any monomial is taken to be of the form b^a'^ Q. If, besides the polyno- 
mials, all entire functions in a, b are considered, then the extended universal 
enveloping algebra of the Heisenberg algebra appears or in other words, 
the extended Heisenberg-Weyl algebra. In the (extended) Heisenberg-Weyl 
algebra one can find the non-trivial embedding of the Heisenberg algebra: 
non-trivial elements obeying the commutation relations (|8|), whose can be 
treated as a certain type of quantum canonical transformations. We say 
that the (extended) Fock space, T is determined if we take the (extended) 
universal enveloping algebra of the Heisenberg algebra and attach to it the 
vacuum state |0 > such that 

a|0 > = . (9) 

It is easy to check that the following statement holds if the operators a, b 
obey (^), then the operators 

= \?a — nb , 
J0 = &a-^, (10) 



^Sometimes this is called the Heisenberg-Weyl algebra 
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J„ — a , 

span the s/2-algebra with the commutation relations: 

[J0,J±]=±J±, [J+,J-] = -2jO, 

where n G C |^. For the reahzation (^) the quadratic Casimir operator is 
equal to 

c^2-^R+,J„-}-« = -^(^ + i) , (11) 

where {,} denotes the anticommutator and is c— number. If n G then 
( |To[ ) possesses a finite-dimensional, irreducible representation in Fock space 
leaving invariant the space of polynomials in b: 

Vn{b) = (1,6,&2,...,6")|0), (12) 

of dimension dim7^„ = (n + 1). The spaces Vn possess a property that 
Vn C Vn+i for each n £ Z+ and form an infinite fiag and 

It is evident that any polynomial in generators J^'" operator preserves the 
flag of V . Such an operator we will call SI2- exactly-solvable operator. 
Take as an example the s/2-exactly-solvable operator of the form 

1)J- = 46a2-4(6-p-i; 



hf{b,a) = AJ^J--U^ + A{p + -)J- = Aba^ - A{b - p - -)a , 



(13) 

where p is a parameter and J^'^ = J^'^ (see (10)). One can demonstrate 
that the eigenfunctions of /ij are the associated Laguerre polynomials of the 

argument b, Ln ^ (&) and their eigenvalues, En = — 4n. 

As the next step we consider two different realization of the Heisenberg 
algebra (|8|) in terms of differential and finite-difference operators. A tra- 
ditional realization of (^) appearing in all text-books is the coordinate- 
momentum representation: 

d_ 

dy 

where the operator b = y stands for the multiplication operator on y in a 
space of functions /(y). In this case the vacuum is a constant and without 
a loss of generality we put |0 > = 1. Recently, a finite-difference analogue 
of dH) has been found IHI: 

a = V+, b = y{l-5V-) , (15) 

where 



— = dy,b = y, (14) 



^For details and discussion see, for example, 
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is the finite-difference operator, 5 is real number and I'±(— (5) = I?=p((5). A 
remarkable property of this realization is that the vacuum remains the same 
for both cases (^4|)-(0) and it can be written as |0 > = 1. 

Substitution of ( |14D into (|l3|) leads to the operator (5) - the algebraic 
form of the Hamiltonian of the harmonic oscillator. Thus, the operator 
( p!3[ ) can be called the algebraic form of the Hamiltonian of the harmonic 
oscillator in the Fock space. It is evident that the procedure of realization 
of the Heisenberg generators a, b by concrete operators (differential, finite- 
difference, discrete) provided that the vacuum remains unchanged leaves any 
polynomial operator in a, b isospectral. Now let us study another 'face' of 
harmonic oscillator by substituting the realization (ll5| ) in (ll3|). Finally, we 
obtain 

hd{y, D±) = ^[y + 6{p + ^)]D+-A{l + ^)yD. . (16) 

Thus, in the realization by finite-difference operators the corresponding spec- 
tral problem can be defined as 

^[y + 5{p + \)\ Hy + '5) - ^ [(1 + |)y + p + \] Hy) + ^(i + ])y'l'(y - ^) 

= E(t>{y) . (17) 

The operator hd{y,D±) is a non-local, three-point, finite-difference oper- 
ator. It is illustrated by Fig.l. 



(l){y-6) c^{y) Hy + d) 



Fig. 1. Graphical representation of the operator (16) 

In general, the function (/){y) in the rhs of (|l^ can be replaced by (^{y + d) 
or (/>(y — S), or by a linear combination of <j){y it 5), (/>(y). It does not change 
the statement that (17) has infinitely-many polynomial eigenfunctions. 

It is easy to check that 

[yg-5a«]nj ^ yin)j _ 

where = y(y — (^) . . . (^y — nd) is a so-called quasi-monomial and I is the 

identity operator. Using this relation one can show that the eigenfunctions 
of (|l^ ) remain polynomials and furthermore the solutions of (17) are given 
by 

n 

L^r-^\y,6)=Y,af^-^VK (18) 

where a^^"*"^^ are the coefficients in the expansion of the Laguerre polyno- 
mials, Ln~^'^\y) = "^^^Qof^^^^. We call these polynomials the modi- 
fied associated Laguerre polynomials. Simultaneously, the eigenvalues of the 
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equation (16) remain equal to (— 4n), n = 0, 1, 2 . . . and they are the same as 
the eigenvalues of the harmonic oscillator problem (1), (6) and (13). Thus, 
one can say that the operator ( |T6[ ) defines a finite- difference form of har- 
monic oscillator Hamiltonian. 

A natural question can be posed about the most general second-order 
linear differential operator, which (i) has infinitely-many polynomial eigen- 
functions and (ii) is isospectral to the harmonic oscillator (1). Following the 
Theorem |]l| one can show that this operator has a form 

hg{y,dy) = 4(AJ0 + BJ-)J- -4jO + 4(p + i)CJ- = 

4{Ay + B)d^-A[y-{p-^)C]dy (19) 

where A, B, C are arbitrary constants and the generators (^) are realized by 
differential operators (0)- However, by a linear change of variable, y — > ay+ 
P the operator (|l9|) is transformed to (5). Thus, without loss of generality 
we can put A = a = 1 and also C = 1. The eigenfunctions of ([l^) remain 

the Laguerre polynomials but of a shifted argument, L„ ^ (y -|- /?). It leads 
to a statement that among the second-order differential operators there exist 
no non-trivial isospectral deformation of the harmonic oscillator potential 
preserving polynomiality of the eigenfunctions. 



The operator (19) can be rewritten in the Fock space formalism by using 

o 

hg{b, a) = 4(6 + B)a^ - 4[6 - {p - i)]a . (20) 



It is evident that the operator (|20| ) is the most general second order poly- 
nomial in a, which is isospectral to (|l^) and also preserves the space of 



polynomials (12). By substitution (|15|) into the operator (20) it becomes 



transformed into a finite-difference operator 

hg{y,D^) = ABDl + ^[y + 5{p + ^)]D^-Ml + ])yD- (21) 

(cf. ([l^)). Thus, in the realization by finite-difference operators the corre- 
sponding spectral problem can be defined as 

4R 4 1 

_^(y + 26) + ^[y-2B + 6{p + -)] 0(y + 5) 

- ^ [(1 + |)y - I + (p + ^)] <P{y) + ^(1 + ^)y<^(y - S) = E(t>{y) , 

(22) 

and it has infinitely-many polynomial eigenfunctions. 

The operator h(i{y,D±) now becomes the four-point finite-difference op- 
erator, see Fig.2 . 
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Fig. 2. Graphical representation of the operator (|2l|) 

It is quite surprising that a hnear shift of variable y in differential operator 
(|l6|) (which gives nothing non-trivial, see discussion above) leads to occur- 
rence of the extra point in their isospectral finite-difference counterpart. 

Now it is time to ask what would happen if in the expressions (13), (20) 
the operators a, b are not the generators of the Heisenberg algebra (^ but 
the generators of the g— deformed Heisenberg algebra 

[a,b]q = ab — qba = 1, (23) 

where g is a parameter. Following the Theorem proved in Q, one can demon- 
strate that within the g— deformed Fock space built on using g— deformed 
Heisenberg algebra (23) there exists the flag of linear spaces of polynomials 
in b, V (see (12)), which is preserved by the operators (13), (20). By a simple 
calculation one can find the eigenvalues of (13), (20) 

= -4{n} , n = 0, 1, . . . , (24) 

where 

{"} = -; — r ' 



is a so-called g— number and {n} — > n, if g ^ 1. If the parameter q the 
spectra of (13), (20) are real. 



The algebra ( |23| ) has a realization in terms of discrete operators (see, for 
example, |2|) 

a = Vq, b = y , (25) 

where 

^ X f{Qy) - f{y) 
^qfiy) = — n — ^ • 
y(i - q) 

This realization has a property that the vacuum remains the same as well 
as for the cases (p^)-(|T5|) and it can be written as |0 > = 1. Now we can 
substitute ( p5|) in (13) the following operator emerges 

^)J- = 4yp2_4(y_p_ 1; 



h,{y, V,) = U'J- - 4J0 + 4ip + -) J- = 4yp2 _ 4(y - p - ^)P, 



(26) 

where the generators = ba, J~ = a have the same functional form as 
( p!o[ ) but obey the q— deformed commutation relation 

[J^'Jli/, = J''J---^J^J- = -J- , 

forming the g— deformed Borel subalgebra b{2)g of the g— deformed algebra 
sli2)g. 
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The operator hq{y,Dq) is a non-local, three-point, discrete, dilatation- 
invariant operator illustrated by Fig. 3. 



Fig. 3. Graphical representation of the operator (16) 



The spectral problem for the operator ( p^ has a form 



yqiq - 1)^ yq{q - 1)2 y{q - 1 

= • (27) 



or, the rhs can be taken as 



or as 



E^iUiqy) ■ (28) 



= E^'^Uiq'y) ■ (29) 

If for the case ( p7| ) the eigenvalues are given by (^ while for (|8|), ( |29| ) 
the eigenvalues are equal to 

E^^^ = -4g"M , n = 0,l,... (30) 

4?) = , n = 0,l,... (31) 

correspondingly, while in the limit g' ^ 1 all three expressions coincide. The 
spectral problems (|27|)-(p9|) can be considered as a possible definition of a 
g— deformed harmonic oscillator. In the literature it is known many other 
definitions of the g— deformed harmonic oscillator. Such a situation reflects 
an existence of an ambiguity appearing when a q— deformation is performed 
and absence of clear criteria, which can remove or reduce this ambiguity. 
For instance, in the literature it is exploited three types of the g— Laguerre 
polynomials (see, for example, an excellent review [Q]), but it is not clear 
why other possible deformations of Laguerre polynomials are not studied. 

Substitution of (25) in (20) gives a slight modification of the expressions 
(26)- (27). Unlike translation-invariant case it does not lead to a change of 
the number of points in the operator (26) as it happened for the operators 
(16) and (21). 



■^For example, any term in non-deformed expression can be modified by multipliers of 
the type and it can be added extra terms with vanishing coefficients in the limit g — > f 
like (1 - g)^ 
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